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Abstract. We construct "stochastic mappings" between power law probabil- 
ity distributions (PD's) and Gaussian ones. To a given vector N, Gaussian 
distributed (respectively Z, exponentially distributed), one can associate a vec- 
tor X, "power law distributed", by multiplying X by a random scalar variable 
a, N = aX. This mapping is "invertible" : one can go via multiplication by 
another random variable b from X to N (resp. from X to Z), i.e., X = bN 
(resp. X = bZ). Note that all the above equalities 

mean "is distributed as". As an application of this stochastic mapping we 
revisit the so-called "zero-th law of thermodynamics problem" that bedevils 
the practitioners of nonextensive thermostatistics. 
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1. Introduction 

Statistical Mechanics' most notorious and renowned probability distribution (PD) 
is that deduced by Gibbs for the canonical ensemble PG], usually referred to as 
the Boltzmann-Gibbs (BG) equilibrium distribution 

/ ■ \ exp (-(3Ej) 

(!) Pg(i) = y , 

with Ei the energy of the microstate labeled by i, (3 — l/fc^T the inverse tem- 
perature (T), ks Boltzmann's constant, and Zbg the partition function. The 
exponential term Fbg = ex P {—(3E) is, of course, called the BG factor. Recently 
Beck and Cohen [Hj have advanced a generalization (called "superstatistics" ) of this 
BG factor, assuming that the inverse temperature (3 is a stochastic variable. They 
effect a multiplicative convolution to obtain a generalized statistical factor Fas i n 
the fashion 

r°° dd 

(2) F GS = / -f ftf) exp {-(3E) = f o F BG , 

Jo P 

where /(/3) satisfies 

/oo 
-oo 
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Note that the o— sign is used to denote the multiplicative convolution between 
two PDs: multiplicative convolution of two PDs fx and fy is the PD of the product 
of the two corresponding random variables X and Y . 

As stated above, (3 is the inverse temperature, but the integration variable may 
also be any convenient intensive parameter. Superstatistics would arise as a "statis- 
tics of statistics" that takes into account fluctuations of such intensive parameters. 
Beck and Cohen also show that, if f((3) is a \ 2 distribution, nonextensive ther- 
mostatistics (NEXT) is obtained, which is of interest because NEXT, or Tsallis' 
thermostatistics, is today a very active field, perhaps a new paradigm for statistical 
mechanics, with applications to several scientific disciplines 0J|SJE1- In working 
in a NEXt framework one has to deal with power-law distributions, which are cer- 
tainly ubiquitous in physics (critical phenomena are just a conspicuous example 
[Jj). Indeed, it is well known that power-law distributions arise quite naturally in 
maximizing Tsallis' information measure (q is a real parameter called the "nonex- 
tensivity index") 



Tsallis' PDs are encountered in analyzing a rather large variety of physical systems 
[El 03 EH ^2 E], which encourages people to continue investigating the nonextensive 
formalism along multiple viewpoints and a multitude of paths. 

In such a spirit, and further pursuing along the road first travelled by Beck and 
Cohen we show here that, for a fixed temperature T (respectively, any adequate 
intensive parameter r) , there exists a mapping between power law PD 's on the one 
hand, and exponential or Gaussian PD 's on the other one . Using this mapping one 
can transform a Tsallis PD into a Gaussian (resp., exponential) PD and vice- versa 
via multiplicative convolution with a chi random variable of variance unity. 
This mapping nitidly reveals, for the first time as far as we know, the (perhaps 
surprising) intimate relation that exists between the conventional Shannon-Gibbs- 
Boltzmann statistics and NEXT. As a first application we will tackle the vexing 
(for NEXT) problem that revolves around thermodynamics zero-th law that was 
first revealed by Raggio and Guerberoff J3| ( see also E3I El EH Ell EH] > a by no 
means exhaustive list, and references therein). 




subject to appropriate constraints. In the case of the canonical distribution there is 
only one constraint, the energy E, i.e., (E) = K (K a constant) and the equilibrium 
canonical distribution writes f(x) = (l/Z q )[l — (1 — q)/3 q E]^ =T , with (3 q and Z q 
standing for the NEXT counterparts of (3 and Zbg above. It is a classical result 
that as q — > 1, Tsallis entropy reduces to Shannon entropy 



(5) 




2. Mathematical framework 



2.1. Case q > 1. Let us consider a random k— variate vector X = {x%, . . . ,Xk) T 
following a Tsallis distribution with parameter q > 1 (or q < 0). Its distribution 
writes 



(G) 
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with n > k and with the notation x+ = max(0, x). As discussed in fx can be 
obtained as the k— variate marginal of a uniform distribution on the sphere in R™. 
Additionally, a stochastic representation for X can be given in terms of Gaussians, 
which is essential for our present purposes, i.e., 

m x = * 

r 

where N = [Ni, . . . , N n ] T is a Gaussian n— variate vector, N is the k— variate vector 
composed of the first k components of N and r = is a chi random 

variables with n degrees of freedom. We note that the first k components of TV are 
shared by the numerator and the denominator. 

Theorem. If a n is i) a chi random variable with n degrees of freedom and ii) 
independent of a vector X distributed as in @), then the product 

(8) a n X ~ Z 

where Z is Gaussian with identity covariance and symbol ~ denotes equality in 
distribution. 

Proof. See the Annex. □ 

Thus a Tsallis system can be " Gausssianized" simply by multiplying each of its 
components by an independent, scalar chi random variable with n degrees of free- 
dom. For statistical mechanics purposes, the variance is a very important quantity 
because it is proportional to the temperature 0. As defined by |[BJ|, the random 
vector X has covariance matrix Rx = h^ki with 1^ the (k x fc)-unit matrix. Thus, 
the random vector Z, as defined by © has covariance matrix (we use the notation 
EX = (X)) 

(9) EZZ T = EalEXX T = nX-I k = I k . 

n 

2.2. Case < q < 1. Let us first underline a useful connexion between the 
instances q < 1 and q > 1: 

Lemma, if X is distributed according to then 

(10) Y= ■ X = 

Vi - X T X 

is distributed as follows: 

(11) fr(Y)<x (1 + Y T Y)'^ 

with m — n — k, and thus is a k— variate Tsallis random vector with parameter 
-£jp2 < q < 1 • A stochastic representation of Y writes 

, ^ X 

(12) Y =— > 

Am 

where N is a k— variate Gaussian vector and Xm is a chi random variable indepen- 
dent on N with m = n — k degrees of freedom. Note moreover that the covariance 
matrix of Y writes 

(13) EYY T = ENN t eAt = — !— 4. 

Xm m-2 
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Now we can state our main theorem: 

Theorem. ifY is distributed according to j 11)) and ifa n is a chi distributed random 
variable with n degrees of freedom independent on Y , then random variable 

(14) Z = °" V 

Vl + Y T Y 

is Gaussian with unit covariance. Moreover, random variable 

(15) b m = °" 

is chi distributed with m degrees of freedom. Note that b m is dependent on Y. 
Proof. See the Annex. □ 

2.3. Extension of the above theorems. The precedent considerations can be 
extended from the Gaussian case to the one of a purely exponential factor. A given 
Tsallis-distributed vector X gets transformed, via multiplicative convolution with 
a chi random variable, into a vector that is distributed according to an exponential 
PD. In other words, from a Gaussian vector G, one can create a Tsallis distributed 
vector X by multiplying by a random scalar variable. This mapping is invertiblc: 
a Gaussian distributed vector can be obtained by multiplying a Tsallis distributed 
one by another properly chosen random variable. 

2.4. Application: the merging of two independent systems. Suppose we 
have two independent systems (in the sense that their mutual interaction is neg- 
ligible) whose states are described by two Tsallis random vectors, independent as 
well, X and Y with, say, q > 1. If we consider the system Z T = [X T ,Y T ~\ we 
immediately realize that it is not Tsallis distributed since its distribution writes 

(16) fz (A, Y) = (1 - X T X)^- 1 (1 - Y T Y) , 

and, consequently, can not be expressed as a function of A T A + Y T Y (except in 
the Gaussian case g = 1) [131 H3J [13 HU [T7| . 

The problem can, however, be circumvented as follows: if we now " pre- multiply" 
system X (resp. Y) with the same chi-distributed random variable a n as defined 
in theorem 12. II then the system 



(17) Z = a n [X T ,Y T f 

is the merging of two Gaussian systems. Moreover, the covariance of Z writes 

(18) EZZ T = I 2k , 

so that a n X and a n Y are both uncorrelated and Gaussian, and thus independent. 
Finally, considering a new chi distributed random variable b m with m degrees of 
freedom and independent of Z, the Tsallis distributed vector 

(19) x TsalHs = f = ^[X T ,Y T ] T , 

Orn Om 

- and NOT i|16[l - turns out to be the "true" representative of the "merging" of two 
nonextensive systems. Such a vector is not characterized by a fixed temperature T 
(or, more generally, by a fixed value of an appropriate intensive system's parameter 
r), but by a superposition of temperatures "centered" at T (resp., superposition 
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of the intensive parameter centered at r) , exactly in the spirit of superstatistics. 
We can then speak of a "stochastic" zero-th law that would apply for nonexten- 
sive thermostatistics: given two independent systems A, B of equal temperature T, 
the pertinent temperature of the associated composite system A + B at equilibrium 
fluctuates around T . 

3. Physical considerations 

For didactic reasons it is convenient at this point to briefly review the well-known 
connection between variance and temperature for a Gaussian PD before probing 
further into the zero-th law. 

3.1. Quadratic classical Hamiltonians. We consider first the "normal modes" 
hamiltonian for N independent particles in three dimensions 

3N 6N 

(20) H = - cUl + c p mP L = J2 c ^ A 

and consider the location of the many-body system in its 6N— dimensional T— space. 
This is represented by a 6N— vector. Obviously, T— position has probabilities pro- 
portional to a function of PH. We consider two instances: conventional BG ther- 
mostatistics and Tsallis' one. 

• For BG statistics, location is represented by a vector Z distributed accord- 
ing to the (1 / Z) exp (— P H) PD. On account of l|20l) then, Z is Gaussian- 
distributed. 

• For Tsallis' statistics, locations is represented by a vector X that is Tsallis- 
distributed. 

If Z = (zi, . . . , z 6 n) t is BG-distributed, then its covariance matrix EZZ T is a 
diagonal matrix. Because of the equipartition theorem pp , the expectation value of 
each member of the diagonal equals ksT/2. Thus, EZZ T — (ksT /2)I§n , and the 
covariance matrix is proportional to the temperature. 

In the case of Tsallis statistics, let the vectors in (f lfif) be distributed according to 

k 

f x (x) cx [l-(l-g)/3$>?] 

i=l 

(21) Mv) « [l-(l-g)/?]T>? L 

1=1 

so that we are merging two system of inverse temperature f3 = 1/fcsT. The repre- 
sentative of the merged system is (I19f) . Note that Z there is a Gaussian vector of 
temperature T, and Xxsaihs a- Tsallis vector that, although not having a definite 
temperature, possesses a temperature distribution "peaked" at T. The zero-th law 
of thermodynamics for Tsallis distributed vectors entails that if we join two nonex- 
tensivc systems of temperature T and wait until equilibrium is re-established, we 
obtain a total system whose temperature "fluctuates" (a la Beck-Cohen) around 
this temperature T. 
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3.2. "Normal modes" or diagonalizable Hamiltonians. Our approach can 
be extended to more general types of Hamiltonians of the type 

6N 

(22) H = Cizf% pi an integer power, 

i=i 

provided the mixing variable a n is modified as follows. 

i 

Theorem. If A is a diagonal matrix whose element Ai_i = a p ^ where a is i) a chi 
random variable with n — k + 2^ i _ 1 ~ degrees of freedom and ii) independent of 
component distributed as in ij^)) . then the product 

(23) AX ~ Z 
where Z is distributed as 

(24) f z (Z) oc exp 

Proof. See the Annex. □ 

Following now an identical line of reasoning as that of Subsection 2.4, we can 
generalize the stochastic zero-th law obtained there to all classical systems described 
by the above type of Hamiltonian. 

4. Conclusions 

We have here generalized the concept of superstatistics devised by Beck and Cohen 
Its typical "convolution path" that involves fluctuations of intensive parame- 
ters allows one to express generalized factors Bqf hi terms of a superposition of 
exponential PDs. We show the the path can be "inverted", at least with reference 
to Tsallis PDs, in the sense that a normal (or an exponential) distribution can also 
be cast as a superposition of Tsallis PD. Stochastic invertible mappings between 
power law distributions and Gaussian (or exponential) ones have been introduced 
and discussed. 

As an application we have shown that the zero-th law of thermodynamics is obeyed 
by Tsallis PDs only in a stochastic sense. In nonextensive thermostatistics, if two 
systems of temperature T are merged, the ensuing equilibrium total distribution 
exhibits fluctuation around that temperature. 

5. Annex: Proofs 
5.1. Proof of theorem H3J ^From 3.471.3], we have 

(25) exp (—cx) — — j—r / exp (— c/u) u~ a ~ 1 (1 — ux)" du 







I» 

Now replace x by X)j=i x 1 ano - u by w _1 in i12ol) to obtain 



( k \ c a f + °° ( k \ 

exp\-c22 x ij = fjaj J exp(-cw)i; Q - 1 ^l-^-^^J dv 



r(o) 



a-1 

dv. 
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Now, as a classical Gamma integral, 

" +0 ° T {a + k/2) 



r°° / x a-l+Aj T(a + k/\ 



so that 



c 



a+k/2 



f* ^ = FT TT Avi exp (~ cu ) ^ 1+ 2 
T (a + fc/2) 

is the distribution of a chi-squared random variable with k + 2a degrees of freedom 
and variance a 2 = I /2c. Thus 

1 V' /2 ( YH-i^\ c Q I> + fc/2) /- +oc M /1\V, V-*?V _1 

eX P " o~2 = ; 7 'k / a X2a+k(v)[-7=) 1-2^ — 



27T(T Z 



1=1 



dv. 

+ 



C r X2a+fe [v) &xju (xi,...,Xk)d,V 



X^5 

where denotes a k— variate Tsallis distribution obtained as the marginal of the 
uniform distribution on the sphere S n 

g{ ) r (a + k/2) {i _ x t x) <*-1 

(2na 2 c)^T(a) 

where 

n — k 

— - — = a «-> n = 2a + k. 

5.2. Proof of theorem 12.21 We know from the lemma that 

X 



Y = 



Vi - x T x 

and we deduce that 

X 



^l + Y T Y 

^From theorem 12 .11 we know that a n X is normal. Thus 

On 



y/1 + Y T Y' 



-.Y 



is normal. Moreover, the distribution of random variable b m = a n /\J\ + Y T Y ca 
be computed as follows 

hm « I (i + r T r)' 1/2 y- exp p 2 (1 + yTy) ) (i + ^ ^ 

so that 6 m is a chi random variable with m degrees of freedom. 
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5.3. Proof of theorem 13.21 Starting form equality (|25|) . we obtain 



expl-c^a^J = Y{a)j 6XP ^ * ( 1 ~ ^ ( ~t ) ) dv 



X j 

,;=i vu* 



r(a) 7 



1=1 u / \ t=l v u 



+ 



with i = 5Z i=1 — ■ Following the same path as in proof of theorem 13.21 we 
deduce that vector with i— th component a Pi Xi, where Xi is distributed as 

(26) = (1 - a*')"" 1 > 

is distributed as vector Z with 



(27) M,)oce X p^f 
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